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Abstract. The dynamics of fluctuating electric field structures in the edge 
of the TJ-II stellarator, that display zonal flow-like traits, is studied. These 
structures have been shown to be global and affect particle transport dynamically 
pQ. In this article we discuss possible drive (Reynolds stress) and damping 
(Neoclassical viscosity, geodesic transfer) mechanisms for the associated E X B 
velocity. We show that: (a) while the observed turbulence-driven forces can 
provide the necessary perpendicular acceleration, a causal relation could not 
be firmly established, possibly because of the locality of the Reynolds stress 
measurements, (b) the calculated neoclassical viscosity and damping times are 
comparable to the observed zonal flow relaxation times, and (c) although an 
accompanying density modulation is observed to be associated to the zonal flow, 
it is not consistent with the excitation of pressure side-bands, like those present in 
geodesic acoustic oscillations, caused by the compression of the E X B flow field. 



PACS numbers: 52.25.Os,52.30.-q,52.55.Hc 
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1. Introduction 

Mass flows are an active research topic in magnetic confinement fusion due to 
their importance for plasma stability and confinement. In recent years, both 
experimental and theoretical efforts have been made to improve understanding of 
the momentum transport mechanisms that determine the observed plasma rotation 
profiles in tokamaks (see, e.g., [2]). In stellarators without a direction of symmetry 
in the magnetic field strength, the neoclassical non-ambipolar fluxes are formally 
dominant in determining the equilibrium radial electric field and plasma rotation. 
Mass flows in such a system experience a magnetic viscous damping in all directions, 
which tends to keep their amplitude low [3]. 

Fluctuating, flux surface-constant zonal flows have a specific interest, for they are 
thought to be instrumental in turbulent transport control and to be a natural state 
of the drift-wave type turbulence. The question: to what extent different magnetic 
configurations can host zonal flow structures? has gathered some attention [4J 03 |6] . 
For non-symmetric systems, Sugama and Watanabe [U found the suggestive result 
that neoclassical optimization could lead to reduced ZF damping and a turbulent 
transport optimization as a 'welcome side-effect'. Previously, this had been found 
experimentally in LHD's inward-shifted configuration [7]. 

Experimental evidence of long range correlations (LRCs) in distant electric 
potential signals was first found in the CHS stellarator [5] with a double Heavy Ion 
Beam Probe system and later in the TJ-II stellarator with a double Langmuir probe 
system [9]. Shortly aferwards other devices reported similar findings (see [10] and 
references therein). 

In this article we continue the anaysis of the spatio-temporal characteristics of 
zonal flow- like floating potential structures initiated in Ref.[l|. In that reference, a 
method to extract the ZF component from a set of LRCed signals was introduced 
and tested. The reconstructed structure showed actual traits of a zonal flow as 
a collective, fluctuating and transport- regulating structure (see next section). Our 
concern in this work is with the dynamics of these structures: their drive and damping 
mechanisms. We investigate the involvement of the Reynolds stress in the generation 
of the zonal-flow structures and the neoclassical collisional damping and geodesic 
transfer as possible damping mechanisms. 

2. Experimental set-up description 

Experiments where carried out in the 4-period, flexible, low shear stellarator TJ-II 
using electron cyclotron heated plasmas. The typical magnetic field strength is ~ 1 
T and the rotational transform ^-(a) = 1.65 at the last closed flux surface (LCFS) 
(typical minor radius a is 20 cm), dropping slowly to a central value of j-(0) = 1.55. 
Electron temperatures range from 1 — 0.8 keV at the centre to an edge value of 100 — 50 
eV. The floating potential measurements presented in this work were obtained with a 
2D array of 5 x 4 tungsten pins, with radial and poloidal separations of 6 and 3 mm, 
respectively (see figure [1} . Simultaneous measurements are taken with another pin in 
a different sector of the machine, having a long parallel connection length with the 2D 
probe (> 100 m).The sampling rate of the floating potential signals is 2 MHz. 
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Figure 1. Pin arrangement in the 2D probe and orientation with respect to the 
magnetic field structure of TJ-II. 



3. Analysis of the dynamics of the zonal-flow like events 

In previous work [I] we introduced a method to extract the spatial and temporal 
characteristics of the zonal flow- like component of the long-distance correlated floating 
potential signals. This extraction is based on the Singular Value Decomposition. Here 
we give a brief explanation of the results of that analysis and refer the interested 
reader to the original reference for details. 

A large matrix M if formed defining its entries as 

Mij = (f>i(tj) , 

where 4>i(tj) is the i-th floating potential evaluated at the j-th time sample. Here 
tj = (j — l)At and (At) -1 = 2x 10 6 Hz. The size of M is thus m x n, where m is 
the number of floating potential signals (20 + 1 in our case) and n is the number of 
samples in the time window under analysis (6 x 10 4 in our case). 
The SVD of M can be written as the sum 



M 



5> Q u Q K)t 



(i) 



where {u a } a =i...m and {'v a } a =i...n are orthonormal vector sets and are called topos 
and chronos respectively. The spatial information is coded in the topos set whereas 
the temporal evolution is contained in the chronos set. The singular values a a are 
non- negative and are sorted in decreasing order, i.e. [a\ > o~i > • • • > 0). Their 
squared values are the contributions of each of the modes a to the total signal energy, 

i.e. E=z ij \Mt j )f^j: T :^ {M) ^ 

In [1] we tentatively identified the zonal flow component as 

* Z f'=^iu 1 (v 1 )t , (2) 

This identification was motivated by the empirical observations that this perfectly 
correlated spatio-temporal structure was: (a) Low- he, (b) broadband spectrum but 
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Figure 2. Long range correlated, zonal flow-like floating potential modulations 
for three TJ-II shots. Each profile is conditionally average to Azf values in the 
different regions of its pdf as shown in the inset. The conditionally averaged 
distant probe potential is shown with squares. The bars in the left figure 
correspond to 1 standard deviation. 

dominated by low (< 10 kHz) frequencies (no clear peak indicative of GAMs was 
observable) (c) explained most of the observed long-range correlation. Consistently, 
we then showed that (d) it was collective, i.e. fluctuations occurred in the 2D probe 
and the distant probe simultaneously and (e) the zonal flow amplitude defined as 
AzF(tj) = {^Yitj) dynamically modulated the outward particle transport as shown 
by all H-alpha monitors around the device. We refer the reader to the cited reference 
for a more detailed explanation of these points. 

3.1. Conditionally averaged E r excursion and dynamical parameters 

In this section we make use of the extraction technique presented above to study 
the temporal and spatial characteristics of the ZF structure. Regarding the spatial 
structure we are interested to know what is the radial profile of the potential 
modulation and also to verify our definition (Eqj2]) showing that the modulations 
are flux-surface collective. This is done in figures [2] and [3] 

Figure [2] shows the spatial shape of the floating potential profile modulations for 
three different discharges. Each radial point in the profiles is obtained as the poloidal 
average of the floating potentials conditionally averaged to Azf values in different 
regions of its probability density function (PDF, shown in the inset). The different 
colours code different deviations from the mean, i.e. most likely, profile shown in 
black. The modulations in the distant probe are plotted with square symbols located 
around the radial position of the peak LRC. In these shots the distant probe was held 
fixed. It displays modulations of the same sign though somewhat lower amplitude as 
those in the 2D probe. 

In figure [3] we show how the floating potential excursions typically (conditionally 
averaged to Azf/o~ > 1) proceed in time. Again, the collective nature of the 
excursions is observed in the similar waveforms of the 2D and distant probe evolutions. 
A temporal asymmetry is noticeable in these structures. In the next sections we 
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Figure 3. Typical evolution of the floating potential excursions shown in Fig[2] 
Each trace correspond to the poloidal average of the floating potential signals at 
each radius. This spatial average is then de-offset and conditionally averaged to 
Azf/& > !)• Similar evolution is seen in the distant probe signals, showing the 
collective character of the excursions. 



obtain the dynamical parameters associated with the time evolution of E x B velocity 
and compare the result with turbulent Reynolds stress and Neoclassical viscosity as 
possible drive and damping mechanisms. 



3.2. Relation to turbulent Reynolds stress forces 

A typical E x B velocity excursion is shown in figure 01a. This is a derived quantity 
from the floating potential evolution shown in Figj3] for the shot TJII#26240 by 
subtracting the traces at p = 0.86 and p = 0.90 to get an estimate of the radial 
electric field at p = 0.88. We model the evolution as the competition of a drive term 
R(t) and a magnetic viscous damping term, i.e. 

^-Ve = R(t) - vV E . (3) 



Dynamics of zonal flow-like structures in the edge of the TJ-II stellarator 



6 



We can get an estimate of the damping time v~ x by fitting an exponential to the Ve 
decay in figure 01a. The clear peak in the typical floating potential evolution (FigJ3]) is 
indicative of a sudden change in the dynamics (e.g. the forcing) of the structure. We 
therefore assume that flow is freely decaying and the forcing R vanishes on average 
(recall this waveform is conditionally averaged) for t > 30/is. The so obtained decay 
times is v~ x = 25/xs, similar to the one obtained in [1] and in biasing relaxation 
experiments [H] . This value is then used to integrate Eq|3] to infer the acceleration 
R(t) assuming the viscosity is constant in time (figlUb). A time-dependent viscosity 
has been derived in |12j which shows that the effective collision frequency involved in 
the viscosity is augmented by the addition of the growth rate of the ZF. Using this 
simplifying hypothesis we ignore this correction, and note that in general this leads 
to an underestimation of the inferred acceleration. 

A natural question is what physical mechanism provides the required force or 
acceleration. The pin distribution of 2D probe allows to compute the turbulent flux 
of poloidal (or more precisely perpendicular) momentum T Vg = VE. r VE,e (he. the 
perpendicular Reynolds Stress (RS)) at the location of the six central pins. We can 
then estimate the RS acceleration, given by the divergence of the flux, as 

d 

dRS = --T- \VE,rVE,e) (4) 

where, the radial derivative is approximated by finite differences and the poloidal 
average (•) is the mean of the poloidally separated r„ e measurements at each radial 
position. Figure SJb shows the conditional average of ans{p = 0.88) over the same 
time windows used to get the Ve evolution shown in |Ha. There is no observable 
causal relation between the locally measured clrs and the Ve excursions. The inset 
in the figure shows the probability density function of a/jg, which nevertheless shows 
that the inferred acceleration fall within the support of the PDF. In other words: we 
observe RS boosts capable of providing the required acceleration, were they global, i.e. 
representative of the flux surface averaged RS, but observe no temporal correlation 
between those boosts and the Ve excursions. The analysis of the two other shots casts 
similar conclusions. 

We emphasise that this conclusion is subject to several experimental 
approximations in the calculation of the Reynold stress acceleration. First, we are 
implicitly taking a poloidaly and toroidaly localised measurement to be representative 
of the flux surface average. Second, we are approximating the fluctuations in the 
plasma potential by those of the probe floating potential, thus neglecting electron 
temperature fluctuations. 



3. 3. Comparison of relaxation time with the Neoclassical collisional damping rate 

We begin this section by introducing a basic model of flow and viscosity. We follow 
Ref.|13j in the notation and choice of Hamada coordinates (V, 9, £). In this coordinates 
the magnetic field is written as 

B = 2TT X '(V)e e + 27ri;'{V)e i , (5) 

where x an( i V' are the poloidal and toroidal magnetic fluxes respectively and prime 
(') denotes derivation with respect to their argument (the volume enclosed by the flux 
surface V). The s-species flows are assumed to be dynamically incompressible and 
tangent to flux surfaces. They are compactly written as 

u s = E s (iP)e e + A s (i>)B , (6) 
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Figure 4. (a)Typical Ex B velocity excursion at p = 0.88 for shot TJII#26240, 
showing a decay time of ~ 25/j.s. (b) Comparison of the inferred forcing (Eq[3{ 
with the measured Reynolds stress acceleration. The probability density function 
of the RS acceleration shows the events capable of providing excursions like the 
one in (a) or larger. 

in terms of two flux constants E and A. The former can be written in terms of 
the s-pressure p s and electric potential </> as E s (tp) = 2tt ( P z ^ + (j)'(ip)j . With this 

definition, and the property eg x B = (valid for any magnetic coordinate system) 
the perpendicular E x B plus diamagnetic flow is recovered form Bx (EqJBJ). The 
parallel part of eg is the Pfirsch-Schlueter flow, as the Hamada poloidal covariant base 
vector satisfies V • eg = and (eg ■ B) = (^o/2tt)It ( = for a zero toroidal current 
stellarator). The A s term is then the so-called 'bootstrap' flow A s = (u s ■ B) / {B 2 ). 

In our extraction of the zonal flow dynamical parameters, namely acceleration 
and relaxation time, we assumed a freely decaying flow and fitted an exponential to 
get a damping rate. Consistently with this interpretation, we would expect the decay 
of E r = —<p'(r) to be determined by the return of neoclassical currents to an ambipolar 
equilibrium only delayed by the ion inertia through the ion polarization current. To 
make this more explicit we start from the momentum balance equation summed over 
species 

Q 

mi—nui + V • IT + V ■ II e + = j x B , (7) 
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where we keep only the ion's inertia (m e <C nij). Projecting this equation on eg and 
taking flux surface averages we get 

m~ (e e • nui) + (e e ■ V ■ Hi) + (e e ■ V ■ II e ) = = , (8) 

where the quasi-neutrality condition V ■ j = has been used to get the last identity. 
The first term on the LHS of Eql8]is the ion polarization current, while the second and 
third terms are recognized as the ion and electron neoclassical non-ambipolar currents. 

In the above discussion electron viscosity is usually neglected, for n e /IT ~ 
^Jm e /mi -C 1 for equal ion and electron temperatures. While this approximation 
may be correct in ion-root plasmas, we need to keep electron viscosity and neoclassical 
flux to model our low density electron-root plasmas. In this parameter region, the ion 
and electron neoclassical fluxes display a similarly strong E r dependence around its 
equilibrium value |14j and therefore both contribute to the return to the neoclassical 
ambipolarity. 

The projection of the neoclassical pressure tensor is approximated in terms of the 
flow components: 

(eg ■ V ■ n s ) = fi s p E s + fj, s b A s , (9) 

where /x* and ix£ are the poloidal and bootstrap viscosities of respectively. We refer the 
reader to [15] for further details. Analytical expressions for the viscosities coefficients 
in the plateau regime (the relevant regime for the edge of TJ-II plasmas) can be found 
in Rcf. |16j . These expressions were found to be in good agreement with the numerical 
solution of the drift-kinetic equation with DKES in LHD. A similar agreement has 
been found in TJ-II [T7]- In this reference it is also shown that the poloidal viscosity 
term in Eq|9] is dominant in the edge of TJ-II. Using this approximation in Eq[8] we 
get 

dE 

m t n (eg ■ eg) = -(jii* + pQ(Ei - E io ) +R . (10) 

Here, EiQ = e ^ i (v'i + p' e ) represents the ambipolar ion poloidal flow obtained 
by zeroing the time-derivative of Eq(5] R represents other possible momentum fluxes 
(i.e. Reynolds' stress) not captured in EqJU] Therefore, our estimate of neoclassical 
damping rate is given by 

nun (eg ■ eg) 

A usual approximation of the inertia term (eg ■ eg) for large aspect-ratio and circular 
cross-section is (eg - eg) w r 2 (l + 2q 2 ). In the case of TJ-II this approximation 
is only correct within a factor of 2. The poloidal viscosity is related to the non- 
ambipolar radial diffusion and was computed with DKES |18j : the mono-energetic 
neoclassical coefficients were calculated following [19] and then convolved with a 
Maxwellian distribution describing the equilibrium plasma. Indeed, it can be seen 
that EqlHl are equal to the equations proposed by Sugama and Nishimura [T!5] for 
poloidal mass flow balance neglecting the parallel and poloidal heat fluxes (i.e. keeping 
only the particle flux terms). More precise calculations are underway. Also note that 
the calculated viscosities correspond to a fixed value of E r , namely the ambipolar 
electric field calculated self-consistently [50]. We therefore assume that the viscosity 
nonlinearity is not too strong, as is to be expected in the low poloidal Mach number 
conditions discussed here. Finally, we must again note that the relatively fast (~ 10/is) 
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Figure 5. Comparison of NC decay times with the observed ZF relaxation 
time. Calculations for three magnetic configurations (labeled with their volume 
and rotational transform at the last closed flux surface) are shown. The actual 
configuration used for the shots in fig[2] is in blue. The radial decay of v v is 
mainly due to the decaying T-profiles. For the for these calculations we took 
Tj (0.9) = 50 eV and T e (0.9) = 120 eV. 



time-scales involved in the rise and decay of the potential structures would require a 
timc-dependent treatment of the viscosity (i.e. without neglecting the explicit time 
derivative of the distribution function in the drift- kinetic equation (T2l |4"]V 

The comparison of the estimated NC decay times with the experimental values 
is shown in Figure El Despite the aforementioned approximations, a reasonable 
(within a factor x2) agreement is found. A more systematic comparison would 
require local temperature and density measurements (to be used as an input for the 
viscosity calculations) and the extension of these studies to magnetic configurations 
with different damping characteristics, like those shown in FigJSJ This is left to future 
work. 

3.4- Density profile modulations and comparison with density compression by the 
E x B flow 

In Ref.pQ we reported on the correlation between the ZF amplitude and the H a 
monitors as a proxy of the outward flux of particles. Positive ZF amplitude 
(corresponding to a positive increment of E r ) was seen to be dynamically correlated 
with decreased levels of H a radiation. This correlation is shown in figure |6]a (inset) 
for the shot no 26971. In this shot the 2D probe was operated with columns 1, 2 and 
4 (see figure [1]) operated in floating potential mode and column 3 set to ion saturation 
current (I sa t) mode. The extraction method outlined above was applied to the floating 
pins together with the distant probe potential to obtain the <pf modulations shown 
in Fig|6]a with estimated radial electric field {—^4>f) arid I sa t modulations shown in 
Fig(5]b and c. The line averaged density is somewhat higher than for the shots shown 
in figure [5] and the average floating potential profile displays more negative values 
approaching the electron to ion root transition. 

Ion saturation current depends mainly on plasma density and electron 
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Figure 6. Floating potential (a), approximated radial electric field (b) and Ion 
saturation current (c) modulations associated to the LRCed ZF-like structure in 
the edge of TJ-II. The inset in (a) shows the anti-correlation between the H a 
radiation and Azf first reported in pp. The inset in (c) shows the values of 
twice the geodesic curvature in the probe region. These were used to compute the 
expected E X B compressional density modulations shown in gray. 



temperature I sa t oc n\fT\., so its modulation is to be interpreted as modulations 
of these quantities with a dominant contribution from the density. We observe a 
steepening (flattening) of the I sa t profile for positive (negative) values of Azf- This is 
consistent with the observed anti-correlation C(Azf, H a ). However we note that, in 
this shot, the floating potential profiles observed for negative Azf amplitudes, exhibit 
large E x B shearing rates around p = 0.85, when and where the density profile 
flattens. 

As shown in figure EJa, the ZF rise times are in the range of a few microseconds. 
Pressure waves transit time is qR/c s ~ 10/xs so it is appropriate to ask whether the 
observed modulations in density can be due to compressional effects. Compressional 
pressure sidebands and their associated parallel pressure gradients have been 
considered a mechanism of energy transfer (so-callled 'geodesic transfer') from the 
zonal flow back to turbulence [211 E2] • To first approximation the density sidebands 
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riity, 9, <j>) are expected to be given by the compression of the background flux-constant 
density n (ip), i.e. 

dni 7F 
— — « -n V • v E v E ■ Vni , 
at 

with the second term in the RHS being a factor n\/no smaller than the first. 
Compression of = Bx ^f — is due to toroidicity V • v^; « — 2v£-Vln_B w — 2v£-k, 
where fc is the field line curvature and the approximations arc good for low-/? plasmas. 
Note that the compression of the mean E x B velocity is compensated by the Pfirsch- 
Schlueter parallel flow, so we only consider the compression of the ZF velocity field. 

The inset in figure Elc shows the values of 2n g = 2 gj^ffi ■ k, which according to 
the above discussion relate to the increments of density by 

^«2«^A* ) (12) 
no ts 

where E^ F = — d ^ dr . The values of 2n g are positive in the outer probe region ~ 0.5 
m _1 go to zero around p = 0.85 and turn slightly negative further inside. Figure |Blc 
shows the comparison of the actual density modulation with the E x B compression 
estimate given by Eq[l2]with At = 20/is (see Figj4j) shown in gray. This estimate casts 
considerably smaller relative modulations (a few percent) compared to the observed 
~ 25%. It is well known (see [23]) that the amplitude of density perturbations in 
geodesic acoustic modes are much smaller than their electric potential counterpart. 
Another prominent disagreement in this comparison is that the location of maximum 
E r modulation and expected density compression coincides with the radial position 
where I sa t hardly varies. 

All the above observations indicate that the observed density modulations are 
not due to the compression of the background density by the E x B flow, but rather 
to radial transport of particles in agreement with our previous findings pQ. 



4. Summary and Conclusions 

In this work we have studied the dynamics of the zonal-flow like structures responsible 
of the long range floating potential correlation observed in the TJ-II stellarator under 
certain conditions [24]. We made use of a previously introduced [T] ZF extraction 
method to discuss possible drive (Reynolds stress) and damping (Neoclassical viscosity, 
geodesic transfer) mechanisms for the associated E x B velocity. We show that: (a) 
while the observed turbulence-driven forces can provide the necessary perpendicular 
acceleration, a causal relation could not be firmly established, possibly because of the 
locality of the Reynolds stress measurements, (b) The calculated neoclassical viscosity 
and damping times are comparable to the observed zonal flow relaxation times. The 
configuration dependence of this decay will be addressed in future work, (c) Although 
an accompanying density modulation is observed to be associated to the zonal flow, 
it is not consistent with the excitation of pressure sidebands, like those present in 
geodesic acoustic oscillations, caused by the compression of the E x B flow field. 
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